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Abstract 

We present spectroscopic observations for 64 radio galaxies having ultra steep 
radio spectra. Twenty-nine objects have redshifts 2; > 2, the largest redshifts 
being almost 4. Our ultra steep spectrum (USS) criterion {a < —1) has proven 
to be the most efficient way of finding distant radio galaxies. We find that 
even among the USS sources, there is a strong statistical correlation between 
the spectral index and redshift. The most distant radio galaxies within the USS 
sample have the steepest radio spectra. 

In our sample there are 3 radio galaxies at z > 3 compared with 26 at 2 < 
z < 3. However, the present data do not allow us to decide whether there is a 
decrease in co-moving source density at the highest redshifts. 

We have analyzed the spectra of the 30 objects with the highest redshifts 
{z > 1.9). For these high redshift radio galaxies, Lya is almost always the 
dominant emission line, with a rest frame equivalent width ranging from ~ 100 
A to more than 1000 A. The equivalent widths of the most important emission 
lines (Lya, CIV, Hell, ClII]) are found to correlate strongly with each other. 
The large rest frame equivalent widths and the correlation between the equivalent 
widths of the emission lines, con&m that photoionization by a central continuum 
source is most likely the dominant ionization mechanism. 

There are significant velocity differences between the various emission lines of 
our high redshift radio galaxies; in particular the Lyct line is shifted with respect 
to the higher ionization lines. Velocity shifts range from 100 to almost 1000 km 
s^^ in some cases. Simulations show that the effects of associated HI absorption 
on the Lya emission line may be responsible for most of these velocity shifts. 
However, other mechanisms such as organized kinematics of the Lya emission 
line gas (e.g. inflow or outflow) and obscuration of the line emission from the far 
side of the radio galaxy may also play a role. 

Key words: Galaxies: active - Galaxies: redshifts 
1 Introduction 

The most distant observable galaxies are radio galaxies. Because these objects can be seen 
to large redshifts and because they are spatially extended, they are unique probes of the 
early Universe and important laboratories for testing models of galaxy formation. The 
most efficient way of obtaining significant samples of radio galaxies at redshifts ^ > 2 is to 
concentrate on the counterparts of radio sources with ultra-steep spectra (USS, 5* ~ z/", a ^ 
— 1.0; e.g. Miley & Chambers 1989; Chambers Sz Miley 1989). During the last 5 years more 
than 60 radio galaxies with z > 2 have been discovered (including those described in this 
paper), mostly using this technique. The most distant USS radio galaxy discovered so far is 
8C1435+633 at z = 4.25 ( |Lacy et al. T9Mi| ). 

High redshift radio galaxies (HZRGs) exhibit a large variety of properties. Remarkable is 
the initially unexpected alignment between the UV/optical emission and the radio structures 



2. Samples 



page 2 



and the common presence of enormous (> 100 kpc) ionized gas halos for which the Lyman a 
emission hne is very strong (e.g. McCarthy et al. 1993). The variety of emission components 
from HZRGs and their large luminosities provide unique diagnostics for studying galaxies in 
the early Universe. 

In order to study the nature of HZRGs and their evolution, it is important to obtain a 
large sample of such objects. After the initial success of the high redshift search technique 
on a sample of 4C USS sources (8 of 33 had z > 2, Miley & Chambers 1989; Chambers & 
Miley 1989), we set about increasing the number of HZRGs by extending the USS search 
method to larger and fainter samples of radio sources. The resultant Leiden intermediate-flux 
compendium of USS sources was the basis for a key programme carried out at the European 
Southern Observatory in Chile between 1990 and 1994. This programme was remarkably 
successful, resulting in the discovery of 22 new radio galaxies with z > 1.9. 

The project consisted of several steps. First, samples of USS sources were compiled having 
a range of finding frequencies and flux densities (see Rottgering et al. 1994. Secondly, high 
resolution (1.5") imaging of USS sources with the VLA was carried out (see Rottgering et 
al. 1994) to obtain good positions so that the USS sources could be identified on deep CCD 
images. Thirdly, a subset of the sources was imaged in R band with 2-m class telescopes 
to obtain reliable identifications (see Rottgering et al. 1996a).) Finally, low-resolution 
spectroscopy was carried out with 4m-class telescopes on a subset of the objects to determine 
their redshifts. Here we present the results of this spectroscopy. We have augmented it 
with previously unpublished data from the original 4C sample in order to carry out some 
initial statistical analyses. The spectra presented here have been the starting point for 
several follow-up programmes to investigate various properties of HZRGs, some of which are 
discussed elsewhere in this thesis. 

The structure of this paper is as follows. In Section 2, we describe the sample selection 
and in Section 3 the observations and reduction of the data. The redshifts and spectra are 
presented in Section 4. We also analyze the properties of the emission lines in this section, 
determining velocity shifts between the emission lines and searching for correlations with 
other properties of the radio galaxies. In Section 5 we discuss our results. We consider the 
efficiency of the USS technique for finding high redshift radio galaxies and briefiy discuss 
the radio luminosity function of USS sources at high redshifts. Furthermore, we investigate 
a possible origin for the observed velocity shifts of Lya with respect to the higher ionization 
lines as being due to associated HI absorption systems. In Section 6 we summarize our 
results and conclusions. 

Throughout this paper we assume a Hubble constant of Hq — 50 km s~^ Mpc~^ and a 
deceleration parameter of qo — 0.5. 

2 Samples 

The sources whose optical spectroscopy we shall discuss here, have been selected either from 
the compendium of USS sources as compiled by Rottgering et al. (1994) or from the 4C 
sample of USS sources as presented by Tielens et al. (1979). 
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2.1 Leiden USS compendium 

The Leiden compendium of USS radio sources (Rottgering et al. 1994) contains sources 
from a variety of radio catalogues, including the 8C at 38 MHz ( Rees 19^01) , the 6C at 151 



MHz ( IHales et al. 1988|) , the Texas Sky Survey at 365 MHz (Pouglas et al. 1980|) , the 



Molonglo Catalogue at 408 MHz ( [Large et al. 1981| ), the 1400 MHz NRAO Sky Survey 



Condon ( [Condon fc Broderick 1985| ; [Condon fc Broderick 1986[) and the 4.85 GHz Sky 



Survey ([Condon et al. 1989[) . The resulting compendium therefore contains sources with a 
range of finding frequencies. The median flux density at 365 MHz is ~ 1 Jy. A large fraction 
of the sources from this compendium has been imaged at -R-band to a level of i? > 24 using 
2-m class telescopes at ESO and La Palma (e.g. Rottgering et al. 1996a) 

Since our main aim was to obtain a sample of z > 2 radio galaxies for subsequent de- 
tailed studies, we selected sources for spectroscopy so as to optimize the chance of finding 
high redshift objects. In carrying out the project, compromises had to be made in the se- 
lection of targets. It was regarded as of paramount importance to increase the number of 
high redshift radio galaxies. In order to optimize the use of the scarce 4m-class telescope 
time towards this goal, additional selection criteria were applied to reduce the number of 
spectroscopic targets. While recognizing that these additional criteria introduce complica- 
tions on subsequent statistical analyses, we regarded it as more important at this stage to 
increase the number of high redshift galaxies whose properties could be studied in detail. 
Two additional radio-based criteria were therefore used in optimizing the spectroscopic tar- 
gets for high redshift objects. First, we excluded the counterparts of diffuse radio sources. 
Diffuse (Fanaroff- Riley Class 1, Fanaroff & Riley 1974) radio morphologies are known to be 
associated with relatively weak (nearby) radio sources. Secondly, we concentrated on the 
counterparts of smaller (< 20") radio sources. Taking account of the well-known correlation 
of radio angular size with redshift for quasars ([Barthel fc Miley 1988[) , only 25% of 2; > 2 
steep spectrum quasars have angular sizes > 10". Hence, even if radio galaxies are system- 
atically larger by a factor of 2, this criterion should exclude relatively few z > 2 objects. 
Finally, a selection on optical magnitude from the CCD imaging was applied, concentrating 
on the objects with faint optical identifications {R > 21) as the best candidates for distant 
radio galaxies. 



2.1.1 4C USS sample 

The 4C USS sample comprised 33 sources from Tielens et al. (1979) that have spectral 
indices a < — 1 between 178 and 5000 MHz. Further details on the sample selection and 
observations can be found in Chambers & Miley (1989); Chambers, Miley and van Breugel 
(1987) and Chambers et al (1996a). Within this sample 8 have z > 2 ( [Chambers et al 




including 4C40.36 at ^ = 2.3 ([Chambers et al. 1988[) , 4C41.17 ([Chambers et al 



rggq ) and 4C48.48 at z = 2.3 ( [Chambers fc Miley 19gg[ ). Here we will present the spectral 
properties for these 8 sources. Their median flux density at 365 MHz is ~ 2 Jy. 
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3 Observations of the Leiden USS sample 

The observations were carried out with CCD sht spectrographs on the 3.6m telescope and 
the NTT (3.5m) of the European Southern Observatory (ESO) at La Silla and on the WHT 
(4.2m) of the Observatorio del Roque de los Muchachos at La Palma. The ESO observations 
were carried out as part of a "Key Programme" . Details of the observational setup are given 
in Table 1. Typical observing times were 1-2 hours per object. 

The emission line regions of distant radio galaxies are typically extended by at least 
several arc seconds and oriented along the radio axis (e.g. McCarthy et al. 1990b). In order 
to optimize the signal to noise for detecting such emission lines, we therefore usually aligned 
the slit in the direction of the radio source and used a relatively large slit width (> 2"). The 
typical emission line widths of distant radio galaxies is ~ 1000 — 1500 km s~^. With the slit 
width and the gratings used (see Table 1), the resolution of the spectra are comparable to 
the expected width of the emission lines. 

Since the targets were too faint to be seen on the TV guiders, the following procedure was 
used for positioning the objects in the slit. The slit was first centered on either a reference 
star from the Space Telescope Guide Star System (CSS) or from a star on the R-band 
CCD frames ( [Rottgering et al. 1996a|) , whose position had been accurately determined with 



respect to the GSS system. A short exposure was made to check that this fiduciary star was 
indeed centered in the slit and the telescope was then offset to the target. Since the relative 
position between the target and the positional reference star is known with an uncertainty 
of less than 0.5", the target was assumed to be in the slit. 



3.1 Reduction and Analysis 

Reduction of the spectra was carried out using the 'Long-slit' package in the IRAF image 
reduction package developed by the US National Optical Astronomy Observatory (NOAO). 
For the data taken on every night, first a bias was constructed averaging 10 "zero second" 
exposures taken at either the beginning or at the end of each night. This bias was subtracted 
from every non-bias frame. The pixel-to-pixel gains were calibrated using fiat fields obtained 
from an internal quartz lamp. Wavelength calibration was carried out by measuring the 
positions on the CCD of known lines from either an He-Ne or a Cu-Ar calibration lamp, 
fitting a polynomial function to these data, and applying the resultant calibration factors. 
We estimate the accuracy of the wavelength calibration to be about 5 % of the resolution, 
i.e. ~ 1 A (see Table 1). 

The sky contribution was removed from the frames, by subtracting a sky spectrum ob- 
tained by fitting a polynomial to the intensities measured along the spatial direction on the 
two-dimensional spectrum not including the spatial rows/columns where the targets were 
positioned. One-dimensional spectra were extracted by averaging in the spatial direction 
over an effective aperture as large as the spatial extent of the brightest emission line. Each 
night one or two standard stars were observed and used for the fiux-calibration. We estimate 
that the fiux-calibration is good to ~ 10%. 
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3.2 Identification procedure 

After all the geometrical and intensity calibrations were applied and the sky background level 
was subtracted, identification of the various emission lines in the spectra was attempted. 

Depending on their redshifts, the emission line spectra of radio galaxies are dominated 
by one or several of the following lines: Lya A 1216, [O II] A 3727, [O III] A5007, Ha A 6563 + 
[Nil] A 6583 ( [Spinrad 1986| ; Psterbrock 1989| ; [McCarthy 1993|) . For the more distant {z > 1) 



radio galaxies other lines that are characteristic of the spectra, but with fainter intensities 
(5 - 10% of Lya) are CIV, Hell, ClII], ClI], [NelV], Mgll and [NeV]. 

The procedure for determining the redshifts consisted of first attempting to identify 
detected emission lines with the dominant lines and thereafter searching for fainter lines at 
the appropriate redshift. 

An additional consideration is the presence of the Lya "break" in the continuum spectrum 
due to the Lya forest or a decrease in the intensity of the stellar continuum bluewards of 
Lya. The presence of such a break at the appropriate wavelength is an additional argument 
for identifying a line with Lya. 

We note that 1243+036 (z=3.6) is the only object in the z > 2 sample for which we 
have detected only one emission line. For a number of reasons we are confident that the 
redshift is correct, the most important being that [OIII] 5007 emission is clearly detected in 
the infrared at the expected wavelength (cf van Ojik et al 1996a). 



3.3 Analysis of the spectra 

After an initial identification of the lines, each spectrum was analyzed and relevant pa- 
rameters of the lines were determined. The first step was fitting the line with a Gaussian 
profile superimposed on a constant continuum. From this, the peak intensity and the width 
(FWHM) of the line were obtained. 

The effect of noise on each spectrum was determined by calculating the RMS intensity 
in a line-free spectrum, i.e. excluding wavelengths within twice the measured FWHM of 
the emission lines or residual sky lines. The RMS intensities in these resulting spectra were 
calculated in bins of about 25 pixels. The values for the rms in each bin were fitted with 
4th order polynomials. These polynomials then give an estimate of the local rms at each 
location in the spectra. 

The uncertainty in the wavelength of the intensity peak of an emission line is a combi- 
nation of the systematic uncertainties in the wavelength calibration, estimated to be 5 % of 
the resolution, and the uncertainty in the measurement of the peak position, estimated to be 
(FWHM /(2s/n)), where s/n is the signal to noise ratio of the observed emission line (e.g. 
Condon 1989). The weighted average of the redshifts of the peaks of the emission lines has 
been taken as the redshift of the radio galaxies. The uncertainties in the redshift determi- 
nation were calculated in two ways. First, it was calculated from a weighted average of the 
uncertainties in the peak positions of the identified lines. Secondly, it was calculated from 
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the differences between the observed redshift of the source and the redshifts of the individual 
lines. Since some of the lines have velocity shifts with respect to other lines (see below), 
these two errors are not always in agreement. The uncertainty in the redshift, given in Table 
2, is the largest of these two values. Especially the Lya emission line may be affected by 
HI absorption so that the true centroid of the Lya emission line is not well determined by a 
Gaussian fit. This may cause a different redshift for the Lyct line than for the other emission 
lines in the spectrum (see Sects. 4.2 and 5.4), thus contributing to the uncertainty in the 
redshift . 

The deconvolved widths (FWHM) of the lines were determined assuming that the lines 
are Gaussian. We note that from observations of the Lya emission at resolutions of a factor 
10 higher than the observations presented here it is clear that the line profiles can be very 
complex ( [van Ojik et al. 1996a| ; [van Ojik et al. 1996 b| ). The deconvolved width presented 



here is therefore only a rough estimate of the true width. The formal uncertainty in the 
deconvolved width (FWHM) is a combination of (i) the uncertainty due to the signal to 
noise ratio of the line, estimated to be (FWHM/(s/n)) (e.g. Condon 1989) and (ii) the 
uncertainty in the resolution, estimated to be 10 %. If the fitted FWHM is less than the 
resolution or the determination of the deconvolved FWHM is less than its error, then the 
error in the deconvolved FWHM is taken as the upper-limit to the deconvolved FWHM. 

The line flux was obtained by summing the intensities of the pixels above the determined 
continuum over a wavelength range of size 4 times the width (FWHM) of the line and centred 
at the peak of the emission line. The error in the flux determination is calculated taking 
into account (i) the signal to noise ratio of the line and (ii) the error in the flux calibration, 
estimated to be 10%. We note that in general the emission region will be significantly more 
extended than the width of the slit and that therefore the fluxes of the lines are a lower limit 
to the total flux of the emission line. 

The rest-frame equivalent width of each line was determined as the ratio of the line 
intensity to the local continuum intensity divided by {1 + z). Since the continuum emission 
from our targets is faint, in a significant number of cases only an upper limit to the continuum 
could be determined, yielding a lower limit to the equivalent width. 

In general the spatial extent of the Lya emission line was observed to be significantly 
larger than the point spread function determined by the seeing. The angular size was deter- 
mined by measuring the distance (in arcseconds) between the most extreme spatial points 
where Lya was detected on the two-dimensional spectrum, i.e. at a limiting surface bright- 
ness in the spectra of typically 10"^*^ erg s~^ cm~^ arcsec"^. This measurement is clearly 
dependent on parameters such as seeing and signal to noise ratio of the line and therefore 
the tabulated values should be taken only as an indication of the size. Furthermore, it is 
clear from narrow band imaging that the morphology of the Lya emission is usually complex 
and non-spherical. Hence, extensions in directions other than the slit direction (along the 
radio axis) will differ from the tabulated values presented here. The higher ionization lines 
were usually not detected with enough signal to noise to decide whether they were spatially 
resolved. They may be similarly extended as the optical continuum emission from the radio 
galaxies (usually 2-3"). 
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4 Results 

4.1 Redshifts 

From our spectroscopic observations we have found redshifts for a total of 64 USS radio 
sources, of which 29 are at redshifts larger than 2. In Fig. 1 the individual spectra of the 
objects with z > 1.9 are presented. The lines that we have identified are indicated in the 
plots of the spectra. We also have indicated with the symbol 0, where the spectra are 
severely affected by badly subtracted sky-lines. In a large fraction of such cases we excluded 
the affected pixels and therefore the residual sky-lines are absent on the plots of the spectra. 
In Table 2 the redshifts of the objects with z > 1.9 plus the relevant parameters of the 
emission lines in the spectra are presented. 

Table 3 lists the sources with redshifts smaller than 1.9. In the case of three objects for 
which only one emission line was detected, the redshift needs confirmation and is shown in 
parentheses. The redshift distribution of our sample of radio galaxies is shown in Fig. 2. 

From a few objects broad emission lines were observed (> 4000 km s'^ FWHM) and 
they are therefore probably quasars. The broad line objects with z > 2 are 1113—178 
and 1436+157. The latter object is a quasar with a possible companion galaxy within the 
extended Lya halo of the quasar, oriented along the radio axis. Because no obvious line 
emission was detected from the companion, we cannot definitely say whether this is an 
associated quasar-galaxy pair. 

Statistical evidence for the preferential occurrence of companion galaxies in the direc- 
tion of the radio axis was provided by the R band images of our sample of radio sources 
QRottgering et al. 1996a| ; [Rottgering et al. 1996b| ). However, line emission was usually 
observed only from one of both objects, as in the case of 0828+193 at z = 2.572 and the 
above mentioned 1436+157 at z = 2.538. Therefore for these individual cases it is unclear 
whether both objects are at the same redshift. Only in 1707+105 aX z = 2.349 is the Lya 
emission observed to have two clear spatial peaks corresponding to the R band positions of 
two galaxies oriented along the radio axis. 

4.2 Emission line properties 

In the following we shall investigate the properties of the emission lines and the presence of 
correlations between the derived emission lines properties of the radio galaxies with redshifts 
larger than 1.9 (all objects from Table 2). We do not include the two objects with broad 
(> 4000 km s~^ FWHM) emission lines in the analysis. Thus, the sample for which we 
search for correlations of the properties consists of 28 HZRGs. 

4.2.1 Lya properties 

Lya is the strongest and most spatially extended emission line in the spectra of our HZRGs. 
The measured spatial extent of Lya along the slit ranges from a few arcseconds to 16 arc- 
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seconds, corresponding to ~ 25 to 140 kpc. In Fig. 3 we show the distribution of Lya 
luminosities. The large spatial extent of Lya along the slit indicates that the Lja extent 
perpendicular to the slit may be larger than the slit width (~ 2.5") . Therefore, it is possible 
that the true Lya luminosities are a factor of ~ 2 higher than measured in our spectra. 

Roughly 80% of the radio galaxies have Lya luminosities larger than lO"^^'^ erg s~^. Given 
the sensitivities of our spectroscopy (see Rottgering 1993), this is unlikely to be a selection 
effect. The Lya luminosities are comparable to the radio luminosities of the radio galaxies 
and also comparable to the luminosities of extended Lya emission observed in high redshift 
radio loud quasars ( jHeckman et al. 1991| ). 

Assuming a volume filling factor of the emission line gas of 10~^ as derived from sulphur 
emission lines of low redshift radio galaxies ( [van Breugel et al. 1985| ; pieckman et al. 1982D , 



the Lya luminosities of the radio galaxies imply masses in ionized gas of a few times 10^ to 
a few times 10^ Mq (see also van Ojik 1995a and 1995b). 



4.2.2 Equivalent widths 

The equivalent widths of the various emission lines (Lya, CIV, Hell, ClII]) are plotted 
against each other in Fig. 4. In the cases where the continuum level could not be measured, 
lower limits to the equivalent widths were determined. In a few cases where an emission line 
was not detected but the continuum could be measured, an upper limit to the equivalent 
width of the line was determined. 

We find that the rest frame equivalent widths of all 4 lines are strongly correlated with 
each other. Spearman rank analysis of the correlations give confidence levels of more than 
99%. The strongest and most obvious correlations are visible between the equivalent widths 
of the higher ionization lines (CIV, Hell, ClII]). The correlations with the Lya equivalent 
width are less significant. 

We note that most of the lower limits imply that the true correlations between the 
equivalent widths might be even stronger than seen in Fig. 4. We found no evidence for 
correlations of the equivalent widths with other properties of the radio galaxies within our 
sample (radio size, luminosity, distortion, spectral index). 



4.2.3 Velocity shifts between emission lines 

In many cases we find that the redshifts of the peak intensities determined using Gaussian 
fits to the individual emission lines in a spectrum differ by more than the combined errors in 
the redshifts of the individual lines. The rms-uncertainties in the wavelength calibration are 
typically less than 2 A, corresponding to ~ 130 km s~^, while the measured velocity shifts 
between the lines are between several hundreds and sometimes almost 1000 km s~^. 

To establish the reality of these shifts we must ensure that there are no systematic 
errors in our wavelength determination. An uncertain wavelength determination may occur 
when an emission line occurs at a wavelength outside the range of the lines emitted by the 
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calibration lamp, so that the wavelengths are determined by extrapolation rather than by 
interpolation from the calibration lines. This is the case for 0448+091, where Lya is at 3692 
A and the detections of CIV, Hell and ClII] are marginal. In the observing session during 
which 0448+091 was observed (WHT January 1991) there were few calibration lines in the 
far blue. The bluest usable calibration line was at 4131 A. Thus, the Lya emission line was 
340 A blueward of the nearest calibration line. This extrapolation may well have caused 
an erroneous determination of the wavelength of Lya and we therefore ignore this galaxy in 
considering the velocity correlations. During the other observing sessions there were no such 
problems and calibration lines were available down to 3888 A. 

Finally, there is the problem is the presence of strong skylines. In the objects 0529—549, 
0828+193 and 4C 26.38 the CIV emission line is affected by the 5577 A skyline. We have 
therefore excluded these three objects in the following analysis of the velocity shifts. 

After also leaving out 1243+036 for which Lya was the only detected emission line in the 
spectrum, the final sample for which we have carried out the velocity shift analysis consists 
of 23 HZRGs. The determined velocity shifts between the various emission lines and their 
uncertaintes are listed in Table 4. In this table negative values indicate that the first line 
is redshifted with respect to the second. Although the uncertainties on most velocity shifts 
are large and the shifts have only a 2-3 a significance, the presence of velocity shifts in the 
sample is statistically significant (see below). 

In Fig. 5 we show the velocity shifts between the emission lines Lya, CIV, Hell and 
Cm] plotted against each other. Several correlations between the velocity differences are 
apparent: strong correlations are present between the measured velocity differences of the 
higher ionization lines with respect to Lya. Objects with larger velocity shifts in one high 
ionization line with respect to Lya tend to show also larger shifts of the other high ionization 
lines with respect to Lya. This is strong evidence that there is a real systematic shift between 
the redshift of the peak intensity of Lya with respect to the other lines. 

From Fig. 5 it can also be seen that in almost all objects CIV is redshifted with respect 
to both Hell and ClII]. The average shift of CIV with respect to Hell is —332 ± 101 km 
s~^ and Civ with respect to Clll] is 437 ± 125 km s~^. There may also be a correlation 
between the velocity differences of C IV with He II and C IV with C III] , although the errors 
on the individual points of the plot are larger than in the case of the first three correlations. 

One object that has a relatively large offset in the plots from the other objects and 
relatively large uncertainties is 1138-262. This radio galaxy is a very peculiar galaxy, both 
from its R band image, where it is an extremely clumpy galaxy (Pentericci et al. 1996| ), and 
from a high resolution radio image ([Carilli et al. 199BD , where it appears as a string of radio 



knots. Furthermore, the high resolution Lya spectrum reveals that it has a very peculiar 
velocity structure ([Pentericci et al. 19961) . Although it has these peculiarities which may be 
the cause of its particular position in the velocity plots, we have no reason to believe that 
our velocity determination has any systematic error as in the case of 0448+091. 

We find from a Spearman rank analysis that the correlations of velocity shifts of the 
higher ionization with respect to Lya are significant at a confidence level of more than 
?o. The apparent (ClV-Hell)-(ClV-ClII]) correlation would have a formal significance 
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in a Spearman rank analysis of 95%, but the errors on the individual points are too large 
to allow a Spearman rank analysis. Furthermore, if the previously mentioned radio galaxy 
1138-262 is omitted, the correlation becomes insignificant (85%). Thus, we conclude that this 
apparent correlation is not statistically significant. No significant correlations were found 
between the other velocity differences. Thus, the three most significant correlations in the 
velocity shifts indicate that there is frequently a velocity shift between the Lya line and the 
higher ionization lines. 

Summarizing, we find that in many cases significant velocity shifts between emission lines 
are present. Velocity shifts of Lya with higher ionization lines are strongly correlated and 
Civ appears to be mostly redshifted with respect to Hell and ClII]. A comparison with 
velocity shifts earlier reported in quasars and radio galaxies ( [Gaskell 1982 ; Wilkes 1984; 



Ihlspey et al. 1989| ; |Corbin 199(i| ; [Uarswell et al. 1991| ; |lVlcCarthy et al. 1990a| ; [I'adhunter 



19|9l|; [Bales et al. 1993b| ) and their possible interpretation will be discussed in a next section. 



5 Discussion 



5.1 USS selection efficiency 

The technique of selecting ultra steep spectrum radio sources, has proven to be very successful 
in finding distant radio galaxies. In total we have obtained spectra of 147 objects, of which 
for 64 (44%) we have been able to determine a redshift. From the remaining objects no clear 
emission lines were observed. There are four possible explanations for the failure to detect 
emission lines in these objects. First, the objects did not have emission lines of sufficient 
strength to be detectable in the given atmospheric conditions and sensitivities. Secondly, 
many of these sources may be at redshifts between ~ 1.2 and ~ 1.6, where there are no 
extremely strong emission lines visible in the optical window. Thirdly, some objects may be 
at z > 5, where Lya is shifted beyond the red of the window to be observable. Fourthly, 
there may have been some unseen errors in the telescope pointing. 

Of the 147 objects observed, 29 (19.7%) are at redshifts larger than 2, which shows that 
the USS selection is indeed an efficient way of finding very distant radio galaxies. Taking 
into account the preselections from VLA imaging of ~ 600 sources and from R band CCD 
imaging of ~ 300 objects, we can say that of the radio sources with spectral index a < — 1, 
radio size < 20" and faint optical identification {R > 21), roughly 1 in 5 objects is at a 
redshift larger than 2. 

To further assess the efficiency of the USS selection technique we examined the relation 
between redshift and spectral index. In Fig. 6 we plot this relation for a total of 108 radio 
galaxies, where the spectral index is that defined between the lowest available frequency 
above 150 MHz and 5 GHz. Apart from our USS sources (selected a < —1), this plot also 
includes radio galaxies from the Molonglo Survey (a < —0.9, McCarthy et al. 1990a, 1990b), 
the Bologna Radio catalogue (a < —0.9, McCarthy et al. 1991; Lilly 1988), 8C ( [Lacy et al.| 
r994a[ ; [Lacy et al. 199^ ) and 6C ([Eales et al. 1993a[) . 



From this plot we see that there is a clear trend for the higher redshift sources to have 
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steeper radio spectra. The median spectral index of all sources is —1.15. At 2; < 1 the median 
spectral index is —1.08, at redshifts between 1 and 2 it is —1.16, at redshifts between 2 
and 3 it is —1.19 and at 2; > 3 the median spectral index is —1.31 A Spearman Rank 
correlation analysis shows that the correlation between spectral index and redshift is present 
at a confidence level of more than 99.9%. If only considering the radio galaxies at redshifts 
larger than 1, this correlation is still significant at the 98% confidence level. 

Particularly, for the steepest spectra the efficiency of the technique in finding the most 
radio distant radio galaxies is remarkable. Of all sources with radio spectra steeper than 
— 1.3 almost 30%±15% are at redshifts larger than 3, compared with only 3%±2% of the 
sources with spectral index between —0.9 and —1.3. 

This continued steepening of the observed radio spectra with redshift can be understood 
as the effects of the "K-correction" at radio wavelengths: There is a well known tendency 
of radio sources to have spectra which steepen at higher frequencies. For objects at higher 
redshifts the receivers sample higher emitted frequencies and steeper spectra. Evidence that 
the intrinsic radio spectra of HZRGs indeed continue to steepen to higher frequencies comes 
from 5 and 8 GHz VLA observations of our sample ( |Carilli et al. 199"6| ). An additional effect 
that could contribute is that the radio emission would have an intrinsically steeper spectrum 
at a fixed rest frame frequency due to increased inverse Compton losses from scattering of 
the microwave background, which is stronger at higher redshifts (Krolik & Chen 1991; see 
also Blumenthal & Miley 1979; Chambers et al. 1990). 

Further it has been suggested that the spectral index might be linked to the total radio 
power ( |Laing fc Peacock 19801 ) possibly due to the density of the environment. Because we 



are observing a fiux limited radio sample, the objects at the highest redshifts may have a 
larger average radio power. However, at 2; > 2 the range in radio power of the sources is 
relatively small. There is no evidence for a correlation between radio power and the spectral 
index of the z > 2 sources (see Carilli et al. 1996). Although a mechanism in which a denser 
environment of a radio source (as might be expected around primeval galaxies in the early 
Universe) may cause the radio lobes to have a steeper spectrum, the most plausible cause 
for the strong relation of spectral index with redshift, remains the steepening of the radio 
spectrum at higher rest frame frequencies with a possible contribution from increased inverse 
Compton losses due to the more intense microwave background at high redshifts. 



5.2 The evolution of the space density of radio galaxies 

Although there appears to be an abrupt decrease in the observed number of radio galaxies at 
redshifts 2 > 3, this does not necessarily imply a real decrease in the space density of radio 
galaxies at those redshifts (e.g. Dunlop & Peacock 1990). A nice way to test for evolution 
beyond a defined redshift zq is to use a modified version of the <V lVmax> test, the "binned" 
<^'I^Lax> (e-g- Osmer & Smith 1980). 

For such a test it is important that the fiux limit of a survey is well established. The 
fiux limits for our samples of USS sources are not well defined for three reasons. First, our 
sample from the Texas Catalogue and the sample from the 4C survey do not have a well 
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determined flux limit. At low flux levels (< 0.5 Jy for the Texas Catalogue and < 2 Jy for 
4C) the surveys are incomplete. Secondly, the actual flux limit for the Texas Catalogue is 
a function of source morphology. Thirdly, the selection on spectral index by combining two 
radio catalogues increases the flux limit at low frequencies for sources with steeper spectra. 

This poorly known flux limit introduces an error in the <V' /Vl^^^> of at least 0.1-0.2 
(see Rottgering 1993). Therefore we conclude that our present sample of USS radio galaxies, 
does not allow us to decide whether there is a decrease in source density at the highest 
redshifts. 



5.3 Ionization of the emission line gas 

Three important ionization mechanisms that have been considered for the emission line 
regions of radio galaxies are shock ionization from the interaction with the radio lobes, 
photoionization by hot stars and photoionization by the active nucleus. In general, the line 
ratios of radio galaxies and the rest frame equivalent width of Lya are not well reproduced by 
shock ionization models, nor by hot stars, and photoionization by a nuclear UV continuum 



appears to be the dominant mechanism of ionization (|Baum & Heckman 1989a; 


Baum & 


Heckman 1989b 


; Baum et al. 1992 


; McCarthy 1993 


; IFerland & Osterbrock 1986 


; Ferland & 


Usterbrock 1987; 


Chambers et al. 199q; |Charlot & l^all 1993 


)• 



The strong correlations between the equivalent width of the emission lines found for our 
sample is further support that the emission lines are all produced by the same ionization 
mechanism and are consistent with central photoionization. 



5.4 Velocity shifts: possible associated absorption systems 

Velocity shifts between different emission lines in the spectrum of quasars occur frequently. 
For the broad lines of quasars, the high ionization lines and Lya are systematically blueshifted 
with respect to the low ionization lines and Ha by as much as 4000 km s~^ ( [Caskell 1982| ; 
Wilkes 1984| ; |Espey et al. 1989| ; |Corbin 1990] ; [Carswell et al. 1991|) . These velocity shifts 



have been attributed to outflow of highly ionized gas from the active nucleus in combination 
with obscuration of the far side of the quasar's emission line region. 

The characteristics of the velocity shifts in our HZRGs are quite different from those of 
quasars. Not only are the shifts generally smaller, but also the shifting of the peak of the 
Lya line with respect to the peak of the higher ionization lines is not observed in quasars, 
where Lya together with the higher ionization lines are shifted with respect tot the lower 
ionization lines and Ha. Furthermore, in quasars the shifts are observed in the broad line 
region, i.e. on a very different scale than the narrow line region of radio galaxies. It is 
therefore plausible that the origin of the velocity shifts between the emission lines in HZRGs 
is different from that in quasars. 

We consider two explanations for the occurrence of the velocity shifts in our sample of 
HZRGs. First, the spatial distribution and the velocity fleld of the high ionization lines and 
the Lya emission lines are not necessarily the same and could lead to the observed velocity 
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shifts. In radio galaxies there have been reports of velocity shifts of up to ~ 1000 km s ^ 
between the emission lines ([McCarthy et al. 1990a| ; rTadhunter iMTt [Bales et al. 1993b|) . 
These have also been proposed to be related to the nuclear activity ( [ladhunter 1991[ ) or to 
be due to inflow or outflow on the near side of the galaxy, with obscuration of the emission 
from the far side of the galaxy ( pi)ales et al. 1995b[ ). 

Secondly, a large fraction (60 %) of HZRG show the presence of strong associated HI 
absorption in spectra of Lya emission at a resolution a magnitude higher than the spectra 
presented in this paper (Rottgering et al. 1995a; van Ojik et al. 1996b). In most cases the 
Lya emission is absorbed at wavelengths blueward of the peak of the Lya emission, resulting 
in a red-shift of the peak of the Lya emission when observed at low resolution. 

Since the observed correlations appear to be mainly due to the Lya emission line shifting 
with respect to the higher ionization lines, this indicates that HI absorption systems may 
play a role in producing the observed velocity shifts. In the next paragraph this will be 
further considered. 

5.4.1 Simulations of Lya velocity shifts caused by HI absorption 

Because of resonant scattering of Lya photons, associated H I absorption systems can produce 
a strong absorption feature in the Lya emission profile that alters the shape of the line. When 
observed at low spectral resolution this may result in a slightly different redshift for the Lja 
line than for the higher ionization lines. An absorber located at a small blueshift (or redshift) 
with respect to Lya will cause a slightly higher (or lower) redshift for the measured peak of 
Lya. To investigate whether the effect of associated HI absorption systems on Lja could 
indeed produce the range of observed velocity shifts, we have simulated the profile which 
would be observed if a Lya emission line with associated H I absorption was observed at low 
spectral resolution. 

We assume that the original emission line profile of Lya is Gaussian. An absorption 
feature is produced in the emission line profile by convolving it with a Voigt profile for a given 
H I column density and offset velocity from the Lya emission peak. The resultant theoretical 
profile of Lya emission with absorption feature is then convolved with the instrumental 
resolution. The peak position of this simulated "observed" profile is determined and gives 
the offset velocity with respect to the peak of the original Gaussian Lya emission line profile. 

We have carried out the simulations for an instrumental spectral resolution of 20A 
(FWHM) which is representative for our observations, and an observing wavelength of 4000 
A for Lya, so that the instrumental resolution is 1500 km s~^. The parameters that influence 
the resultant Lya velocity shift are the offset velocity of the HI absorption system, the HI 
column density, the Doppler parameter (b) of the absorber, and the width of the original 
Gaussian Lya emission line. From the high resolution observations of 0943—242 ( [Rottgerm^ 
3t al. 1995[) it was found that the HI absorption system has a column density of ~ 10^^ 
cm~^ with a Doppler parameter of 55 km s~^, while the velocity width of the original Lya 
emission profile was estimated to be 1575 km s~^. The velocity width of Lya is usually in 
the range of 1000-1500 km s"^ (FWHM) (e.g. McCarthy 1993; see also Section 3 and van 
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Ojik et al. 1995b). Therefore, we shall assume similar parameters of the Lja emission and 
H I absorption systems for our simulations. 

Column densities range from 10^^ to 10^^ cm~^, Doppler parameters of the absorption 
systems range from 25-75 km s~^ and the offset velocity of the absorber is varied from to 
1500 km s~^ We have performed the simulations for input Gaussian emission lines of widths 
1800, 1500 and 1000 km s'^ (FWHM). 

Figure 7 shows the resultant velocity shifts for the observed Lya profile as a function 
of absorber offset. The solid lines are for column densities 10^^, 10^^ and 10^° cm~^ with 
6 = 50 km s^^. For each column density the lower and upper dashed/dotted lines indicate 
the shifts for b = 25 and b = 75 km s~^. The different plots are for Lya velocity widths of 
1800, 1500 and 1000 km s~^ (FWHM). Figure 7 also shows the amount of flux that remains 
after HI absorption as a function of absorber offset, for HI column densities and Lya widths 
the same as above. The solid lines are for 6 = 50 km s~^, while lower and upper dotted lines 
now correspond to 6 = 75 and b = 25 km s~^. Figure 8 indicates the maximum shift of the 
Lya peak that can be produced by H I absorption systems of different column densities and 
Doppler parameters. 

We find from these simulations that, depending on the velocity width of Lya emission 
line, H I absorption systems can produce velocity shifts in the observed peak of Lya emission 
of 100-600 km s~^ for HI column densities of 10^^-10^" cm~^. Such velocity shifts are 
comparable with the velocity shifts observed in our sample of radio galaxies. Thus, the 
velocity shifts between Lya and the higher ionization lines could be explained if the majority 
of HZRGs have strong (N(HI) ~ 10^^-10^° cm~^) associated absorption. 

However, the presence of associated HI absorption systems cannot be the cause for all 
observed velocity shifts between the various emission lines. Of the objects with significant 



velocity shifts, there are at least five with known strong H I absorption in the Lya line ( |van 
O jik et al. 1996bD , but also four that do not have strong associated HI absorption. These 
last objects actually have some of the largest observed significant velocity shifts between the 
emission lines. Furthermore, the origin of the observed systematic shift of C IV with respect 
to Hell and ClII] is unclear. Possibly, different velocity profiles cause the observed shifts 
between these emission lines if they originate from different regions in the galaxies. We note 
that Civ is also sensitive to absorption. However, for an HI absorption with high column 
density (> 10^^ cm~^), the associated CIV absorption column density is usually a factor 
10^-10^ lower (e.g. Bergeron & Stasihska 1988) and will not have a very strong effect on the 
position of the peak of the line. 

Thus, HI absorption systems cannot be the only cause of velocity shifts. Additional 
mechanisms such as outflow or inflow combined with obscuration, as proposed for quasars, 
or radio source induced displacements must be important. 



5.5 Continuum alignment of the z > 1.9 Leiden compendium 



It is well known that the major axes of the optical continuum and the radio emission of radio 
galaxies at z > 0.8 are aligned, while low redshift radio galaxies have a weak tendency to 
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have the optical minor axis ahgned with the radio axis ( [Chambers et ah 1987| ; [McCarthy 



et ah 1987). Interpretations of the "ahgnment effect" include jet-induced star formation 



and scattering of anisotropically emitted nuclear continuum ( [Begelman fc CiofH 1989[ ; [Ree; 
19|8|; [di Serego Alighieri et al.T989| ; rTadhunter et al. 1992[) . 



The R band continuum of many of the z > 0.8 radio galaxies is contaminated by strong 
emission lines extended in the direction of the radio axis. Such strong extended emission 
lines might cause the alignment effect, which is often measured at the lowest flux levels of 
the extended emission. It has therefore been suggested that the observed alignment does 
not involve the optical continuum (e.g. Meisenheimer et al. 1994). At redshifts of ~ 1, the 
strong [O II] A3727 emission line is in R band. This emission line is usually spatially extended 
in the direction of the radio axis (see McCarthy et al. 1991b). At redshifts between 1.9 and 
3, however, there are no strong emission lines in R band. The only significant emission line 
present in R band at these redshifts might be CIII]A1909. From the measured equivalent 
widths measured for our sample of HZRGs we deduce that C III] contributes less than 10% to 
the R band flux and often less than 5%. Although C III] may be slightly spatially resolved, 
it is very unlikely that it would cause a strong alignment. 

For 18 of the z > 1.9 objects of the Leiden USS compendium we have the radio position 
angle from 20cm VLA maps ( K.ottgering et al. 1994] ) and an R band image from which 



we could reliably determine the optical position angle ( [Rottgering et al. 1996aD . Figure 



9 shows the difference between the optical continuum axis and the radio axis for these 18 
objects, divided into two groups with less than 5% line contamination in R band and with 
line contamination between 5% and 10%. The alignment is similarly strong for both subsets, 
indicating that the minor contamination by line emission from C III] in R band does not play 
a role. This plot provides further evidence that the previously inferred alignment between 
the optical continuum of HZRGs and the radio structure is real. 



6 Radio source sizes 



The observed differences between steep spectrum quasars and radio galaxies have been ex- 
plained as merely due to whether the radio source axis is more closely aligned with the line 
of sight (quasars) or in the plane of the sky (galaxies) (see Antonucci 1993 and references 
therein). Clearly this model predicts that the projected linear sizes of quasars should on 
average be smaller than that of radio galaxies. Barthel (1989) found that steep-spectrum 
quasars with 0.5 < z < 1 have sizes that are systematically smaller than radio galaxies with 
0.5 < z < 1 by a factor 2.2, indicating a cone angle of ~ 45° along which the quasar can 
be seen. Gopal-Krishna and Kulkarni (1992) have obtained a similar results for a sample of 
radio sources with 0.1 < z < 2. At low redshift, it seems that there is a relative shortage 



of observed steep spectrum quasars ([Kapahi 1990[) , indicating that the simple unification 



scheme with a fixed cone angle breaks down at low z. 

We have compared the radio source size distribution for the z > 2 USS radio galaxies 
presented here with the radio source sizes for the z > 2 quasars from the quasar sample 
studied by Barthel and Miley ([Barthel fc Miley 1988[) (see Fig. 10). In both cases the 
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linear sizes have been calculated using a Hubble constant of Hq — 75 km s~^ Mpc~^ and a 
decelaration parameter of Qo — 0.5. The mean angular radio size of the quasars is 18 ±6 and 
of the USS radio galaxies 40 ± 7.5. The ratio between the sizes is 2.2, the same factor as 
has been found for the lower redshift objects. Within the context of the model of Barthel, 
this would imply that the cone angle of ~ 45° that had been found for the sources with 
0.5 < z < 1 also holds at high redshift. 

However, there are two selection effects. First, as we have pointed out previously, there 
is an angular size bias in the USS sample in favour of small radio angular sizes. Secondly, 
the radio galaxies all have have been selected to have ultra-steep spectral indices, whereas 
the quasars form a proper flux limited sample with a median spectral index of —0.8. We 
show the projected linear size of the radio galaxies and quasars as a function of spectral 
index in Fig. 11. This shows evidence that the steepest-spectrum quasars tend indeed to be 
the largest. Restricting the sample to USS quasars, the 4 remaining sources are too few for 
a meaning full statistical comparison with the USS radio galaxies. 

We therefore conclude that the difference in the source size distribution for the z > 2 
USS radio galaxies and radio galaxies is merely tentative. A proper comparison between the 
sizes of radio galaxies and quasars should include data covering a similar range of spectral 
indices for each species. 

7 Summary and conclusions 

In this paper we have presented spectroscopic results for galaxies associated with 64 ultra 
steep spectrum radio sources. Of these radio galaxies, 29 arc at z > 2 and the three highest 
redshifts are 3.6, 3.8 and 3.8. Our ultra steep spectrum (USS) criterion (a < —1) has proven 
to be the most efficient way of finding distant radio galaxies. We find that even among the 
USS sources, there is a strong correlation between the spectral index and redshift. The most 
distant radio galaxies within the USS sample have the steepest radio spectra. 

In our sample there are 3 radio galaxies ai z > 3 and 26 at 2 < z < 3. However, the 
present sample of USS radio galaxies does not allow us to decide whether there is a decrease 
in co-moving source density at the highest redshifts, because of the uncertain flux limit of 
the sample of radio sources. 

Wc have analyzed the spectra of the 30 objects with the highest redshifts {z > 1.9). For 
these high redshift radio galaxies, Lja is almost always the dominant emission line, with 
a rest frame equivalent width ranging from ~ 100 A to more than 1000 A. The equivalent 
widths of the most important emission lines (Lya, CIV, Hell, Clll]) are found to correlate 
strongly with each other. The large rest frame equivalent widths and the correlation be- 
tween the equivalent widths of the emission lines, confirm that photoionization by a central 
continuum source is most likely the dominant ionization mechanism. 

There are significant velocity differences between the various emission lines of our high 
redshift radio galaxies; mainly the Lya line is shifted with respect to the higher ionization 
fines. Velocity shifts range from 100 to almost 1000 km s~^ in some cases. Simulations show 
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that the effects of associated H I absorption on the Lya emission hne may be responsible for 
most of these velocity shifts. However, in at least a few objects, other mechanisms such as 
organized kinematics of the Lya emission line gas (e.g. inflow or outflow) and obscuration 
of the line emission from the far side of the galaxy must play a role. 

The sample of objects discussed here is well suited for follow-up observations. Various 
detailed studies of their properties are presented in other papers or are underway, both with 
ground-based telescopes and with the HST. 
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Fig. 1. One dimensional spectra of 2; > 1.9 USS sources 
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Fig. 2. Redshift distribution of the USS radio galaxies in our sample whose redshifts could 
be determined 
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Fig. 4. The rest frame equivalent widths of the emission hues Lya, CiV, Hell and CIII] 
plotted against each other. Lower limits are indicated in the plot as upward and rightward 
pointing arrows, upper limits as downward and leftward pointing arrows 
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Fig. 5. The velocity shifts between the emission hnes Lya, CIV, Hell and Clll] plotted 
against each other 
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Fig. 6. Radio spectral index plotted against redshift for 108 radio galaxies. Filled dots 
are from our USS sample. The spectral indices are determined between the lowest available 
frequency above 150 MHz and 5 GHz. The open circles are also from our USS sample 
with redshifts obtained at Lick Observatory by McCarthy & van Breugel (1994, private 
communication). The crosses (x) are from the Molonglo sample of McCarthy (1990a, 1990b), 
the pluses (+) are Bologna sources (McCarthy 1991, Lilly 1988), the triangles are from the 
8C survey (Lacy et al. 1994, Lacy 1992) and the diamond is 6C 1232+39 (Bales & Rawlings 
1993). 



8. Figures 



page 33 



500 7 




.00 \ 



200 400 600 800 1000 1200 1400 
Offset absorber (km s"') 




200 400 600 800 1000 1200 1400 
Offset absorber (km s"^) 




FWHM = 1800 km s" 



200 400 600 800 1000 1200 1400 
Offset absorber (km s"') 



1.00 f 



0.90 - 



0.80 




0.60 - 



0.50 - 
0.40 - 



FWHM = 1500 km s"' _ 



200 400 600 800 1000 1200 1400 
Offset absorber (km s"') 



1.00 f 




FWHM = 1000 km s"' 



200 400 600 800 1000 1200 1400 
Offset absorber (km s"^) 



Fig. 7. Left) Velocity shifts of observed Lya peaks produced by HI absorption systems as 
a function of absorber offset velocity. Plots are for different N(HI) as indicated (in cm~^). 
Solid lines are for b — 50 km s~^, dashed/dotted lines are lower and upper values calculated 
for b = 25 and 75 km s~^. Right) Flux ratios of Lya before and after absorption with 
indicated HI column density. Solid lines for 6 = 50 km s~^, lower and upper dashed/dotted 
lines for b = 75 and 25 km s~^. For more details see text 



8. Figures 



page 34 



1000 




14 15 16 17 18 19 20 21 

HI column density (cm"^) 

Fig. 8. Maximum velocity shifts of the peak of Lya as a function of HI absorption column 
density. The solid lines are for the indicated original widths of Lya emission (in units km 
s~i FWHM) and b — bO km s~^. Lower and upper dashed/dotted lines correspond to 6 = 25 
and b = 75 km s~^ 
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Fig. 9. Difference between the orientation of the optical R band and the radio axes for 
18 USS sources with 1.9 < z < 3. At these redshifts there are no strong emission hues 
with wavelength within the filter bandpass. The alignment is similarly strong for objects 
with hue contaminations in R band of 5%-10% as for those objects with less than 5% line 
contamination 
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Fig. 10. Comparison of the projected linear sizes for z > 2 radio sources. Left shows the 
steep spectrum quasars from the sample of Barthel and Miley. Right shows the distribution 
for our USS radio galaxies 
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Fig. 11. The linear size as a function of spectral index for the z > 2 USS radio galaxies 
(left) and quasars from the Barthel and Miley sample (right) 
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9 Tables 



Table 1. Observational setup 



Dates 


Sep 90 


Nov 90 


Mar 91 


Nov 91 


Jan 92 






May 91 




Apr 92 












Nov 92 




Telescope 


NTT 


ESQ 3.6 


NTT 


ESQ 3.6 


WHT 


Instrument 


EF0SC2 


EFOSC 


EMMI 


EFOSC 


ISIS 


CCD 


Thompson^ 17 


RCA # 8 


Thompson 


TEK # 26 


EEV 4 


CCD size 


10242 


1024 X 640 


1100 X 1040 


512 X 512 


1280 X 1180 


Binning 


2x4 


2x2 




1x2 


3x2 


Binned Pixel 


0.7 


0.67" 


0.44" 


0.61" 


0.74" 


Slit 


2.2" 


2.5" 


2" 


2.5" 


2.5" 


Grism(s) / gratings 


UV 300 / 


B300 


Grism 3 


B300 


R158B 




B300 










Spectral ranges [10^ x A] 


3.6 - 5.4 / 


3.6-7 


4.1 - 7.9 


3.8 - 6.9 


3.6 - 6.9 




4.7-6.8 










Wavelength resolution 


25 A 


24 A 


18 A 


24 A 


20 A 
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Table 2. Redshifts and emission line properties for z > 1.9 objects 



Name z Line Peak Flux FWHM Eq. width Lya 

Wavelength 10"^^ Extent 









(A) 




(erg s 


^cm ^) 


(km s' 


/ 


(A) 




(") 


0200+015 


2.229+ 0.002 


Lya 


3927 


+ 2 


17.4 


+ 2.2 


< 681 




> 256 




6 






CIV 


5002 


3 


4.2 


0.5 


1658 


624 


63 


13 








Hell 


5297 


2 


3.2 


0.4 


890 


506 


36 


5 








CIIIJ 


6158 


2 


A n 

4.0 


0.5 


1423 


487 


58 


9 




0211-122 


2.336+ 0.002 


Lya 


4053 


+ 1 


5.7 


+ 0.6 


1698+ 


635 


49 + 6 


10 






NV 


4139 


2 


4.1 


0.5 


1812 


666 


35 


4 








SiIV/O 


4674 


6 


0.6 


0.2 


< 736 




5 


2 








CIV 


5168 


1 


5.6 


0.6 


1225 


472 


48 


5 








TT TT 

Hell 


5472 


1 


O 1 

3.1 


f\ A 

0.4 


898 


453 


33 


4 








CTTTl 




1 

± 


9 9 


n 3 




41 1 


27 


3 




0214+183 


n 1 on 1 r\ nno 

2.130+ 0.003 


CIV 


4852 


+ 1 


3.0 


+ 0.3 


1142+ 


516 


81+ 12 








Hell 


5125 


4 


1.8 


0.2 


2016 


586 


50 


8 








CTTTl 


5Q71 





1 8 

-L.O 


n 9 


801 

OU -L 


493 


38 


5 




0355—037 


ci 1 HO 1 n nm 

2.153+ 0.001 


Lya 


3836 


+ 2 


11.2 


+ 1.8 


< 667 




> 263 




8 






CIV 


4883 


4 


2.7 


0.6 


859 


690 


26 


7 








Hell 


5170 


4 


3.7 


0.7 


1514 


687 


62 


18 








PTTTl 


fini X 

uuxo 


Q 
O 


9 




^ 47n 




33 


9 




0417—181 


2.773+ 0.003 


Lya 


4594 


+ 2 


3.0 


+ 0.3 


2303+ 


590 


17 


+ 1 


3 






CIV 


o Aci 

5842 


3 


1.2 


0.2 


1047 


510 


8 


1 








TT TT 

Hell 


6185 


2 


0.5 


0.1 


^ A on 

< 430 




3 







0448+091 


2.037+ 0.002 


Lya 


3692 


+ 1 


12.2 


+ 1.4 


< 558 




> 321 




13 






CIV 


4717 


10 


1.2 


0.4 


1542 


885 


> 101 










Hell 


4986 


4 


1.4 


0.4 


688 


562 


> 40 










/^tttI 

CIIIJ 


5789 


7 


n ^ 

2.7 


0.6 


1709 


930 


> 60 






0529-549 


2.575+ 0.002 


Lya 


4345 


+ 1 


7.4 


+ 0.8 


1944+ 


591 


39 


+ 4 


5 






CIV 


5539 


3 


U.4 


0.1 


< 530 




2 











Hell 


5866 


3 


0.6 


0.1 


< 524 




4 











CIIIJ 


6833 


10 


1.8 


0.3 


2322 1047 


11 


1 




0748+134 


2.419+ 0.004 


Lya 


4161 


+ 1 


6.3 


+ 0.8 


1311+ 503 


> 522 




5 






CIV 


5292 


3 


1.8 


n o 

0.3 


1022 


526 


77 


29 








Hell 


5598 


5 


1.5 


0.3 


1549 


513 


> 209 










CIIIJ 


6515 


6 


1.4 


0.3 


1293 


671 


> 62 






0828+193 


2.572+ 0.002 


Lya 


4342 


+ 1 


13.3 


+ 1.3 


1166+ 548 


> 2730 




7 






CIV 


5534 


1 


1.9 


0.2 


< 435 




43 


5 








Hell 


5862 


2 


1.9 


0.2 


1576 


445 


41 


5 








CIIIJ 


6813 


2 


2.0 


0.2 


1261 


378 


39 


4 




0857+036 


2.814+ 0.003 


Lya 


4640 


+ 1 


2.6 


+ 0.3 


1095+ 441 


113+ 24 


2 






CIV 


5908 


3 


1.0 


0.1 


1650 


426 


78 


17 








Hell 


6247 


2 


0.7 


0.1 


956 


391 


47 


12 
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Name z Line Peak Flux FWHM Eq. width Lya 

Wavelength 10"^^ Extent 
(A) (erg s-icm-2) (km s-^) (A) (") 



DQAS— 242 


2 Q2S+ n 002 




4772 


+ 






CTV 

KJL V 


\J\J 1 \J 


9 






Hell 


6430 


1 






CIII] 


7478 


8 


1113-178 


2.239+ 0.003 


Lya 


3953 


+ 6 






CIV 


5042 


62 






Hell 


5310 


3 






CIII] 


6180 


2 


1138-262 


2.156+ 0.004 


Lya 


3837 


+ 2 






CIV 


4900 


6 






Hell 


5166 


6 






CIII] 


6017 


24 


1243+036 


3.570+ 0.001 


Lya 


5556 


+ 1 


1357+007 


2.673+ 0.002 


Lya 


4467 


+ 1 






CIV 


5682 


5 


1410—001 


9 QRQ-I- 001 




4089 


-1- 1 






CIV 


5209 


1 






HcII 


5516 


2 






CIII] 


6413 


3 




9 00"^ 






in o 






CIV 


5474 


44 






Hell 


5799 


9 






CIII] 


6730 


42 


1545-234 


2.755+ 0.002 


Lya 


4567 


+ 1 






CIV 


5810 


4 






Hell 


6158 


4 


1558-003 


2.527+ 0.002 


Lya 


4287 


+ 1 






CIV 


5468 


2 






HcII 


5785 


5 






CIII] 


6739 


4 


1707+105 


2.349+ 0.005 


Lya 


4076 


+ 1 






Hell 


5480 


3 






CIII] 


6384 


22 


2202+128 


2.706+ 0.001 


Lya 


4506 


+ 1 






CIV 


5736 


5 






HcII 


6064 


14 


2251-089 


1.986+ 0.002 


CIV 


4628 


+ 2 






Hell 


4895 


2 






CIII] 


5696 


2 



20.1 




2 




> 653 






3 Q 




4 


1 61 Q 348 


1 1 3 

J- ±0 


37 




2.7 




0.3 


874 241 


83 


29 




2.3 




0.4 


1397 1053 


> 91 






6.4 


+ 


0.7 


5021+ 764 


51 


+ 3 




1.7 




0.4 


10196 1530 


8 


1 




0.7 




0.2 


< 559 


3 







2.8 




0.3 


1708 459 


16 


1 




13.9 


+ 


1.6 


2603+ 695 


> 832 






0.8 




0.2 


1982 728 


6 


1 




1.3 




0.2 


2787 1005 


10 


1 




1.3 




0.3 


4285 1341 


12 


2 




23.5 


+ 


2.4 


2454+ 251 


351 


+ 13 


7 


9.1 


+ 


1.0 


2188+ 571 


45 


+ 4 


3 


1.7 




0.3 


957 653 


13 


2 




32 




3 9 


19««_l_ 4QC 


356 


+ 57 


10 


5.2 




0.6 


1330 402 


78 


11 




3.6 




0.4 


1617 430 


48 


7 




3.3 




0.4 


1838 444 


73 


13 




42 


_1_ 


4 7 


07011 700 


48 


+ 5 


« 


17.0 




1.9 


9174 1397 


19 


2 




6.0 




0.8 


2952 998 


8 


1 




9.4 




1.7 


4531 1956 


16 


2 




6.7 


+ 


0.7 


1377+ 534 


> 900 




5 


1.1 




0.2 


1221 562 


94 


40 




1.8 




0.3 


1729 560 


> 86 






14.9 


+ 


1.5 


1529+ 556 


1180+ 510 


6 


2.7 




0.3 


1965 501 


95 


15 




1.7 




0.2 


2468 613 


82 


17 




1.2 




0.2 


1389 525 


73 


21 




4.3 


+ 


0.5 


1466+ 501 


242 


+ 80 


16 


0.9 




0.2 


972 475 


70 


24 




0.2 




0.1 


1692 940 


> 7 






7.7 


+ 


0.8 


1167+ 532 


> 1375 




3 


1.9 




0.4 


1464 693 


> 225 






2.4 




0.5 


3023 2071 


> 309 






3.3 


+ 


0.4 


985+ 578 


> 148 






1.3 




0.2 


< 549 


66 


19 




1.5 




0.2 


1016 498 


84 


20 
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Name z Line Peak Flux FWHM Eq. width Lya 

Wavelength 10"^*^ Extent 
(A) (erg s-icm-2) (km s-^) (A) (") 



4C23.56 2.483± 0.003 CIV 

Hell 
CIII] 

4C24.28 2.879± 0.006 Lya 

NV 
CIV 

4C26.38 2.608± 0.001 CIV 

Hell 
CIII] 

4C28.58 2.891± 0.004 CIV 

Hell 
CIII] 

4C40.36 2.265± 0.003 CIV 

Hell 
CIII] 

4C41.17 3.792± 0.001 Lya 

CIV 

4C48.48 2.343± 0.003 CIV 

Hell 
CIII] 

4C60.07 3.788± 0.004 Lya 

CIV 



5388 


± 10 


5.1 




0.6 


5700 


10 


1.3 




0.3 


6651 


4 


1.5 




0.3 


4714 


± 3 


7.3 




1.0 


4835 


12 


6.9 




1.1 


6028 


24 


1.7 




0.4 


5587 


± 2 


8.9 




1.1 


5919 


4 


5.7 




0.8 


6888 


5 


2.4 




0.5 


6033 


± 10 


0.3 




0.2 


6399 


14 


1.6 




0.4 


7424 


5 


1.8 




0.4 


5048 


± 4 


6.2 




0.8 


5356 


3 


5.6 




0.7 


6235 


5 


5.9 




0.6 


5827 


± 1 


22.0 


± 


2.3 


7422 


14 


1.0 




0.5 


5183 


± 3 


6.1 




0.9 


5485 


4 


3.7 




0.6 


6376 


3 


2.8 




0.4 


5831 


± 9 


10.1 




1.3 


7412 


6 


2.7 




0.8 



3695± 1020 149± 38 



1636 


896 


31 8 


734 


626 


34 8 


< 762 




72± 12 


3562 


2149 


70 19 


3284 


2610 


20 6 


< 560 




> 156 


1367 


734 


116 46 


< 603 




> 45 


< 615 




> 1 


2186 


1007 


9 2 


569 


520 


12 2 


2017 


± 782 


54 ± 9 


1640 


718 


54 9 


2370 


668 


67 10 


597 


± 514 


159± 31 


< 685 




> 4 


< 723 




56± 13 


< 766 




44 11 


< 531 




26 4 


2875 


± 937 


153± 71 


< 740 




> 90 
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Table 3. Redshifts of objects from the Leiden USS compendium with z < 1.9 



Name 


z 


Session 


Comments 


0008+172 


1.390 


nov90 




0036+205 


1.370 


aug90 




0050+204 


1.297 


aug90 




0054+090 


1.301 


nov91 




0203-478 


0.836 


nov93 




0241+348 


1.215 


aug90 




0255+114 


(0.447) 


jan92 




0256+324 


(1.657) 


dec89 




0324—228 


1.898 


jan91 




0338-259 


0.440 


nov93 




0410-198 


(0.793) 


nov91 




0429-267 


1.26 


nov91 




/~\ A A ^7 -I A 

0447-164 


1.814 


nov90 




0548-150 


0.650 


nov93 




0729-722 


0.738 


nov91 




0819+672 


0.760 


jan91 




0850+140 


1.106 


jan92 


Quasar 


1035+208 


0.696 


may91 




1043-216 


1.060 


mar 91 




1104+089 


1.384 


jan92 




1228-166 


1.127 


jan92 


Quasar 


1238-273 


1.253 


may91 




1449-004 


1.466 


may91 




1502+039 


1.652 


may 91 




1548-111 


0.959 


apr92 




1635+396 


1.600 


jul91 




1657+003 


1.110 


may 91 




1725+167 


1.508 


may91 




2000-091 


1.238 


may 91 




2214+134 


0.374 


aug90 




2224-273 


1.678 


sep90 




2226+162 


1.519 


nov91 




2245+181 


0.830 


aug90 




2339+269 


0.882 


aug90 
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Table 4. Measured velocity shifts between emission lines 



Velocity shifts (km s ^) 



name 


CIV-Lya 


Hell-Lya 


CIII]-Lya 


CIV-HcII 


CIV-CIII 


Hell-ClII 


0200+015 


-185 + 285 


-83 + 220 


324 + 226 


-101 + 257 


-509 + 263 


-408 + 190 


0211-122 


-528 ± 156 


-315 + 160 


-53 + 156 


-213 + 131 


-474 + 126 


-261 + 132 


0214+183 








-867 + 265 


-623 + 149 


245 + 263 


0355-037 


-19 + 325 


221 + 326 


206 + 233 


-241 + 386 


-226 + 312 


14 + 313 


0417-181 


298 + 218 


532 + 185 




-234 + 190 






0748+134 


295 + 210 


756 + 326 


811 + 315 


-461 + 350 


-516 + 340 


-55 + 421 


0857+036 


-80 ± 177 


503 + 165 




-583 + 207 






0943-242 


-62 ± 148 


279 + 96 


551 + 361 


-341 + 159 


-614 + 382 


-273 + 365 


1138-262 


-968 + 433 


465 + 448 


317 + 1216 


-1429 + 562 


-1281 + 1262 


148 + 1268 


1357+007 


242 + 295 












1410-001 


-197 + 115 


-22 + 152 


330 + 193 


-174 + 158 


-527 + 198 


-353 + 222 


1545-234 


227 ± 237 


91 + 250 




135 + 314 






1558-003 


-599 + 165 


-194 + 276 


-431 + 231 


-403 + 297 


-168 + 257 


236 + 339 


1707+105 




984 + 195 


734 + 1063 






251 + 1071 


2202+128 


41 + 294 


690 + 737 




-648 + 783 






2251-089 








-482 + 210 


-582 + 204 


-100 + 212 


4C23.56 








-443 + 781 


311 + 612 


754 + 576 


4C24.28 


-1336 + 1235 












4C28.58 








362 + 861 


-626 + 563 


-990 + 733 


4C40.36 








456 + 355 


457 + 383 


1 + 333 


4C41.17 


-165 + 573 












4C48.48 








-302 + 293 


-743 + 233 


-441 + 271 


4C60.07 


464 + 564 













